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E865 at the Brookhaven National Laboratory AGS collected about 70,000 K^ 3 events to 
measure the K^ 3 branching ratio relative to the observed K + — > 7r + 7r°, K + — » 7r°/i + i/, and 
K + — > 7r + 7r°7r° decays. The 7r° in all the decays was detected using the e + e~ pair from 7r° — > 
e + e~7 decay and no photons were required. Using the Particle Data Group branching ratios^ 
for the normalization decays we obtain BR(K~^,S) = (5.13 ± 0.02 atat ±0.09 a!/s ±0.04 norm )%, 
where K^ 3 ^ includes the effect of virtual and real photons. This result is w 2.3 standard 
deviations higher than the current Particle Data Group value. Implications for the Vus element 
of the CKM matrix, and the matrix's unitarity are discussed. 



The experimentally determined Cabibbo-Kobayashi-Maskawa (CKM) matrix describes quark 
mixing in the standard model framework. Any deviation from the matrix's unitarity would un- 
dermine the validity of the standard model. One unitarity condition involves the first row 
elements: 

\Vud\ 2 + \Vus\ 2 + \Vub\ 2 = 1-5 (1) 

where a non-zero value of 5 indicates a deviation from unitarity. The V u d element is obtained 
from nuclear and neutron decays. V u b, from the semileptonic decays of B mesons^, is too small 
to affect Eq. (^Q). The V us element can be determined either from hyperon, K — > iriiv(K U 3 ) or 
from K — > Tiev (K e ^) decays. However, K e ^ decays provide a smaller theoretical uncertainty 1 ^'. 
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The most precise value of V u d obtained from the nuclear superallowed Fermi beta decays leads 
to 5 = (3.2 ± 1.4) • 10~ 3 3 , a 2.3d deviation from unitarity. 

Both experimental and theoretical efforts to improv e the d etermination of V u d continue. 
Theoretical contributions to V us were reevaluated recently^^^, but there has been little new 
experimental input on the branching ratio. Since the |V^| 2 and ] V^, s | 2 uncertainties are 
comparable, a high statistics measurement of the branching ratio (BR) with good control 
of systematic errors is useful. 

The bare (without QED corrections) Kj 3 decay rate! ' ' ' ' is: 



dT(K+) = C(t)\V us \ 2 \f + (0)\ 2 [l + A+— ^ (2) 

where t = (P K - P n ) 2 , C(t) is a known k inem atic function, and /+(0) i s the vector form factor 
value at t = 0, determined theoretically Two recent experiments^*^ give A+ (the form 
factor slope) measurements consistent with each other and with previous measurements. An 
omitted negligible term contributing to Eq. |2] contains the form factor /_, and is proportional 
to Mg/M 2 . With current constraints on the hypothetical tensor and scalar form factors'^, their 
possible contribution to the K^ 3 decay rate are also neglibible. 

E865^-^ searched for the lepton flavor violating decay K + — » ir + n + e~. The detector (Fig. 
^) resided in a 6 GeV/c positive beam*^. For the K^ 3 running^*! the intensity was reduced 
by a factor of 10, to 10 7 kaons, 2 x 10 s protons, and 2 x 10 8 7r per 2.8 s spill. The beam was 
intentionally debunched at extraction to remove rf structure at the experiment. The first dipole 
magnet separated particles by charge, while the second magnet together with four multiwire 
proportional chambers (MWPCs: P1-P4) formed the spectrometer. The particle identification 
used the threshold multichannel Cerenkov detectors (CI and C2, each separated into left and 
right volumes, for four independent counters) filled with gaseous methane (Cerenkov threshold 
7t ~ 30 and electron detection efficiency e e ~ 0.98 ^ 3 ), an electromagnetic calorimeter^! anc l a 
muon detector (not used for the measurement). The D and A scintillator hodoscopes gave 
left/right and crude vertical position. 

The 7T° from the kaon decays was detected through the e + e~ from the tt° — > e + e~7 decay, 
with the 7 detected in some cases. To eliminate the uncertainty (2.7%) of the tt° — > e + e~7 BR, 
and to reduce systematic uncertainty, we used the other three major decay modes with a ir° 
in the final state (K + — > 7r + 7r (A' 7 ^ 2 ), K^ 3 , K + — > 7r + 7r 7r (i ; r^3)) for the normalization sample 
(Kdal). 

The K^~ 3 data was collected in a one- week dedicated run in 1998, with special on-line trigger 
logic. 

The Kdal and data were collected by the e + e~ trigger, which was designed to detect 
e + e~ pairs and required at least one D-counter scintillator slat on each (left and right) side 
of the detector and signals from each of the four Cerenkov counters. The Cerenkov efficiency 
trigger required only 3 out of 4 Cerenkov counters (no D-counter requirement). The TAU trigger, 
requiring only two D-counter scintillator hits (one left, and one right), collected events for the 
K + — * 7r + 7r + 7r~ (K T ) sample, to study of the detector unbiased by Cerenkov requirements. 
About 50 x 10 6 triggers were accumulated, ~ 37 x 10 6 in the e + e~ trigger. About 75% of e + e~ 
triggers included accidental tracks, often a jj, from high momentum K — > jj,u or ir — > jiu decays 
partially satisfying the Cerenkov requirement. 

Off-line reconstruction used the spectrometer only. The Cerenkov and D counter efficien- 
cies were obtained from the Cerenkov efficiency triggers. The redundancy of the MWPCs (4 
planes/chamber) and track reconstruction was used to extract MWPC efficiencies. The absence 
of the electromagnetic calorimeter from the trigger allowed its efficiency determination. Each 
efficiency was measured over its relevant phase space. 
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Figure 1: Plan view of the E865 detector with a sim- 
ulated K + — > -K°e + u decay followed by 7r° — » e + e~7. 



Figure 2: Distributions of X and Y positions of the 
first positive track (not e + from the n° decay) for the 
selected K^ 3 and Kdal samples. X and Y positions 
are measured at the end of the second pair of the 
Cerenkov counters (C2). Histograms represent Monte 
Carlo; points with errors represent data. 



Relevant kaon decay chains^^were simulated with GEANtGH (including decays of secondary 
pions and muons). For K+, A+ = 0.0278 ± 0.0019W was used. The radiative corrections to the 
K^ 3 decay phase-space density^ were used. The K+ 3l ( inner bremsstrahlung) decays outside 
the K^ 3 Dalitz plot boundary were explicitly simulated ^1. For ir° — > e + e~j decay, radiative 
corrections were taken into account according to Measured efficiencies were applied^} 

and accidental detector hits (from reconstructed K T events) were added. About 10% of both 
the K^ 3 and Kdal samples had extra reconstructed tracks. 

Selection criteria, common to K^ 3 and Kdal, included requirements for a good quality three 
track vertex in the decay volume (no requirement for exactly three reconstructed tracks was 
applied), for the three tracks to cross the active parts of the detector, for the low (M ee < 0.05 
GeV) mass e + e~ pair to be identified in the Cerenkov counters, and for the second positive 
track to have less than 3.4 GeV/c momentum. The momentum cut rejects events where /i + 
or 7r + from Kdal decays is above Cerenkov threshold and can be identified as e + . A geometric 
Cerenkov ambiguity cut rejected events (27%, 15%, 25%, and 35% for KJ 3 , K^ 2 , K^ 3 , and 
K^ 3 , respectively) where the Cerenkov counter response could not be unambiguously assigned 
to separate tracks^Il 

The K^ 3 sample was then selected by requiring the second positive track to be identified as 
e + in 2 of the 3 electron detectors: CI, C2, or the calorimeter, each with e e ~ 98%. Events 
entering the Kdal sample had no response i n at least one of the two Cerenkov counters. These 
criteria minimized systematic uncertainties but resulted in a small overlap, ~ 3% of the 
K sample and ~ 0.3% of the Kdal which was acco unte d for in the BR calculation. The K^ 2 
acceptance is 1.2%. The K^ 3 acceptance « 0.7%^j somewhat lower because of the lower 
average e + momentum in the K^ 3 decay. The acceptance can be approximately understood by 
assuming a factor of three loss for each charged particle, 30 % for the Cerenkov ambiguity, and 
approximately a factor of 2 for other cuts. Final acceptances for the three modes in the Kdal 
sample differed by < 4% taking into account that either of the 7r°s from K^ 3 can decay into 
e + e~7- The final K^ 3 and Kdal samples were 71204 and 558186, respectively. Figure El shows 
some relevant spatial distributions. 

Contamination of the K^ 3 sample by other K + decays occurred when 7r + or /i + from Kdal 



decays were misidentified as e + , or as a result of ir 
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Figure 3: Energy deposited in the calorimeter by 
the second positive track from the selected K^ 3 sam- 
ple (e + which is not from the low mass e + e~ pair). 
No calorimeter information was used for the e + iden- 
tification. Markers with errors represent data; the 
histogram is simulation. 
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Figure 4: Reconstructed momentum of the e + from 
the low mass e + e~ pair from the selected K* 3 and 
Kdal samples. Histograms represent Monte Carlo; 
points with errors represent data. Plots on the right 
show the bin by bin Monte Carlo to data ratio. 



secondary particle decays was estimated to be at the level of 0.1%. About 8% of final state 
pions decayed into muons inside the spectrometer. The careful MWPC simulation gave good 
agreement of reconstructed track % 2 and vertex distributions between data and Monte Carlo 
calculation. No tight track x 2 cu ts were applied, and the systematic uncertainties estimated 
by variation of the vertex cuts were included in the final result. The check of BR( K T /Kda\), 
described below, also tests the final state tt and \x decays. 

Total contamination of the Ke3 sample was estimated to be (2.49±0.05 s t a t ±0.32 S3/S )%, with 
the systematic uncertainty caused by the simulation accuracy of the CI and C2 response to 7r + 
and fi + . Contamination due to overlapping events was (0.25 ± 0.07)% and (0.12 ±0.05)% of the 
Kdal and K^, respectively. Figure |3] shows the energy distribution in the calorimeter from the 
e + in the sample. The contamination is manifest in the minimum ionization spike at 250 
MeV. The small excess of data in the spike agrees with our contamination uncertainty estimate. 

The final K^ 3 sample included «30% of events with a fully reconstructed 7r°. We used the 
7T° information as a consistency check. Not requiring tt° in our main analysis minimized the 
uncertainty arising from photon detection and reconstruction in the calorimeter, but increased 
vulnerability to contamination from upstream decays and photon conversion. Upstream de- 
cays whose photon produced pairs before the decay volume (evacuated to about 10~ 8 nuclear 
interaction length) were suppressed by requiring the three track vertex to be more than 2 m 
downstream of the decay volume entrance. In addition, the results obtained from the two inde- 
pendent samples, one with and one without the ir° reconstructed, did not show a statistically 
significant discrepancy. 

The Kg S statistical precision is 0.4%. The systematic error estimate, summarized in Table 
^ was determined from the BR stability under variation of reconstruction procedure, selec tion 
criteria, assumed detector efficiencies, and subdivision of both K^ 3 and Kdal samples ^1 No 
significant correlations between any of the different systematic uncertainties were observed. 

The two largest contributions to the systematic error come from the discrepancies^^ between 
data and Monte Carlo in the momentum (Fig. |1J) and spatial distributions. These errors 
were determined by dividing and Kdal events into two roughly equal samples, using the 
relevant parameters, and observing the variation of the result^! The errors were found to be 
uncorrelated. The sensitivity of the vertical spatial discrepancy to the MWPC alignment and of 
the momentum discrepancy to the spectrometer parameters is indicative of their possible origins 



Source of systematic error 


Error esti- 




mate 


Magnetic field uncertainty 


0.3% 


Vertex finding and quality cut 


0.6% 


Vertex position cut 


0.4% 


Cerenkov Ambiguity Cut 


0.3% 


M ee cut 


0.2% 


Detector Aperture 


0.2% 


(7r/^t) + identification 


0.04% 


MWPC efficiencies 


0.2% 


D counter efficiencies 


0.15% 


Cerenkov efficiencies 


0.3% 


Contamination of the selected samples 


0.3% 


Removal of extra tracks 


0.2% 


Vertical spatial/angle distributions discrepancy 


0.8% 


e + /e~ momentum distributions discrepancy 


1.3% 


K^ 3 trigger efficiency 


0.1% 


Uncertainty in the K^ 3 form factor slope 


0.1% 


Total error 


1.8% 



Table 1: Systematic uncertainty 
sources and estimates of their respec- 
tive contributions to the final result's 
uncertainty. The total systematic er- 
ror is the sum (in quadrature) of the 
individual contributions. 



The Z-vertex position is also sensitive to the magnetic field, but has a smaller systematic 
error contribution as determined from both upstream and downstream cuts in Z. 

As an additional consistency check, we estimated the K T /Kda\ BR. T he re sult was (1.01 ± 
0.02) x the Particle Data Group (PDG) ratio^, (the theoretical prediction^! was used for the 
7T° — > e + e~7 decay rate). The 2% error was dominated by the uncertainty in the prescale factor 
of the TAU trigger. A second consistency check compared the K^ 3 BR from 1998 and 1997 data. 
The 1997 K^ 3 data used a trigger that required calorimeter hits, and A and D-counters. That 
trigger neither allowed measurement of these detector efficiencies nor of the trigger efficiency. 
While we did not use the 1997 data for our final result, the 1997 K^ 3 branching ratio was 
statistically consistent (within one sigma) with that from 1998. This agreement is important 
since the momentum spectrum discrepancy between data and Monte Carlo in the 1997 data 
is qualitatively different from 1998 L!j2 '. A preliminary reconstruction version was used for the 
1997 data, without the final magnetic field and detector alignment. This bolsters our intuition 
that the discrepancies in decay product momenta and spatial distributions, which dominate 
the systematic uncertainties, reflect our imperfect knowledge of the magnetic field and detector 
positions but do not bias our result beyond our estimated systematic errors. 

We estimated the form factor slope A + from both 1998 and 1997 K^ 3 data^l We obtained: 
A+ = 0.0324±0.0044 stat for 1998, and A+ = 0.0290±0.0044 stat for the 1997 data, both consistent 
with the current PDG fit. 

After contamination subtraction E9 our result is B i?(!T+, v) / \B R(K^ 2 + K^ 3 + K^ 3 )) = 
0.1962 ± 0.0008 s4ai ± 0.0035 sys , where includes all QED contribut ions (loops and inner 

bremsstrahlung) . As noted above, the 7r° was detected using the e + e~ pair from ir° — ► e + e~7 
and no photons were required. 

Using currentJ2Kdal BR's we infer BR(K+ 3{i) ) = (5.13±0.02 stat ±0.09^ s ±0.04 norm )%, where 
the normalization error was determined by the PDG estimate of the Kdal BR uncertainties. This 
result does not include the correction due to the correlation of the PDG kaon decay ratios, since 
it was estimated to be small compared to the systematic error. The PDG fit to the previous 
K + decay experiments yields BR(K + — > i:®e + v) = (4.87 ± 0.06)% D « 2.3 standard deviations 



lower than our result. 

Radiative corrections for decays inside the K^ 3 Dalitz plot boundary were estimated to be 
— 1.3% using the procedure of Ref.^1 decays outside the Dalitz plot boundary gave +0.5%. 
Thus, the total radiative correction was -0.8% resulting in the bare BR(K+) = (5.17±0.02 stot ± 



0.09 sys ± 0.04 norm )%. 



Using the PDG value for Gp, the short-distance enhancement factor Sew(M„, Mz) = 1.0232 
ECU, and our result for the bare K+ rate we obtain \V us f + (0)\ = 0.2243 ± 0.0022 rate ±0.0007 A 



which gives \V US \ = 0.2272 ± 0.0023 rate ± 0.0007 A+ ± 0.0018 /+(0) if /+(0) = 0.9874 ± 0.0084 ■ 

With this value of V us and V u d from superallowed nuclear Fermi beta decays^, 5 = —0.0003 ± 
0.0016. 

This result is consistent with CKM unitarity, and is i n ag reement with recent theoretical 
determinations of V us from hyperon decays, 0.2250 ± 0.0027^^, and from lattice calculations of 
pseudoscalar decay constants, 0.2236 ± 0.0030^^. However, it increases the discrepancy with 



the value of V us = 0.2176 ± 0.0025 fro m K®^ decay if extracted under conventional theoretical 
assumptions about symmetry brealnng ''. The KLOE K e % measurements for neutral and 
charged kaons currently in progress ^ should help to clarify the experimental situation. 
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